Geometry and topology of knotted ring polymers in an array of obstacles 
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We study knotted polymers in equilibrium with an array of obstacles which models confinement in 
a gel or immersion in a melt. We find a crossover in both the geometrical and the topological behavior 
of the polymer. When the polymers' radius of gyration, Rq, and that of the region containing the 
knot, PiG.k, are small compared to the distance b between the obstacles, the knot is weakly localised 
and Rq scales as in a good solvent with an amplitude that depends on knot type. In an intermediate 
regime where Ra > b > RG,k, the geometry of the polymer becomes branched. When Ra,k exceeds 
6, the knot delocalises and becomes also branched. In this regime, Rq is independent of knot type. 
We discuss the implications of this behavior for gel electrophoresis experiments on knotted DNA in 
weak fields. 
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Circular DNA, as found in bacteria and phages, is of- 
ten knotted Such knots are produced by the action 
of enzymes known as topoisomerases or by random cy- 
clization inside the head of a phage Q. Investigation of 
the types of knots produced can give important insights 
in the working of these processes. The knot type of a 
single DNA can be determined from electron or atomic 
force microscopy. When large numbers of knotted rings 
have to be analysed one uses gel electrophoresis [U] . One 
expects that for a fixed number of basepairs (bp), a DNA 
with a more complex knot (as measured for example by 
the minimal crossing number n c Q) will have a smaller 
radius and hence, when the applied electric field is weak, 
a higher mobility. This expectation is verified experimen- 
tally [f| but despite extensive work Q the precise relation 
between mobility and knot complexity is not fully under- 
stood. For example, the dependence on the length of the 
DNA, as measured by the number of basepairs, has not 
been clarified yet. 

The behavior of melts of polymers of various archi- 
tecture (linear, branched, ring) is also of high interest 
0, [H • In [§] , the geometrical and dynamical properties 
of an unknotted ring polymer in an array of obstacles were 
studied. The array was used as a model for a gel but the 
results were also considered to be relevant for melts of 
ring polymers. Indeed, in a recent experiment (loj . the 
relaxation spectrum predicted in Q was observed in a 
melt of ring polystyrenes. However, since these polymers 
were synthesized near theta conditions, they can contain 
knots This is thus another situation in which topo- 



logical and geometrical constraints may simultaneously 
affect the behavior of a polymer. 



Finally, in the crowded environment of a cell [12j , the 
average size of a macromolcculc is often larger than the 
average distance between macromolcculcs, leading, e.g., 
to important effects on kinetics and geometry of proteins 
[l3j ]. Since some proteins are also knotted 14 1 it is of 



interest to understand how their topological properties 
are modified in crowded in vivo situations. 

Here we study the behavior of a knotted ring poly- 
mer in a regular, cubic, array of obstacles with lattice 
constant b. We remark that besides being a model of 
a gel or melt, such regular arrays can also be realised 
experimentally using hydrogels [1 51 ] and in microfluidics 
Is 12 1 j where they have been used to separate DNA. 



In this Letter, we clarify how the excluded volume con- 
straints exerted by the presence of an array of obstacles 
affect the geometrical behaviors of knotted ring polymers 
in equilibrium. The novel insight concerning knot local- 
ization [l8| and topology dependence of the size of the 
rings will be valuable for the interpretation of gel elec- 
trophoresis of knotted DNA in weak fields. 

Knotted polymer rings are modelled as N edges self- 
avoiding polygons [lH on the cubic lattice. Each edge 
corresponds to one persistence length. In the case of 
DNA this length equals approximately 50 nm or 150 bp. 
For an agarose gel, in which most of the experiments with 
(knotted) DNA are performed, the size of the holes is of 
the order of 200 to 500 nm [j| • In our model, this corre- 
sponds to a 6-value from 4 to 10. In the microlitographic 
arrays of ref. [lq . the effective pore size of l//m cor- 
responds to b Ri 20. Vertices occupied by obstacles, and 
the edges connecting them, are not available for the poly- 
gons. This constraint will decrease the overall entropy of 
the polygons with respect to the free space case affecting 
both their overall conformational equilibrium properties 
and the typical size of the knotted portion. 

The equilibrium properties of the polygons are stud- 
ied by Monte Carlo simulations based on BFACF moves 
[20I ]. This algorithm is known to be ergodic 21 1 within 
the class of fixed knot type and we believe it to be ergodic 
also for polygons confined within an array of obstacles as 
long as b is bigger than 2 lattice spacings. The BFACF 
algorithm works in a grand canonical ensemble where 



2 



configurations are weighted with a step fugacity K [19j . 
The average number of monomers (N) will diverge when 
the fugacity approaches a critical value K c , which is re- 
lated to the entropy per monomer s = —ks In K c [l9j |. 
To increase the Monte Carlo efficiency we implement the 
BFACF algorithm on to a multiple Markov chain sam- 
pling scheme [22[ using several values of K below the crit- 
ical value K c . The values of K c depend strongly on the 
confining geometry and, for each value of b an exploratory 
simulation has to be done in order to estimate K c (b). In 
free space, K c = 0.213538, and from our simulations we 
determined K c (b = 10) = 0.2227, K c (b = 7) = 0.2318 
and K c (b = 5) = 0.2515. As expected K c (b) increases as 
b decreases since for smaller values of the lattice constant 
b, obstacles are more dense in space, the confinement of 
the polymers is stronger and hence the entropy s is lower. 

We next determined the average mean squared radius 
of gyration (R 2 ) as a function of (N) for various b- values. 
We expect a power law relation between these two quan- 
tities, (R 2 ) ~ (N) 2w . In three-dimensional free space, v 
has the value uj — 0.589. We expect that the polymer be- 
gins to 'feel' the presence of the obstacles when its typical 
size Rq = {R 2 } 1 / 2 becomes of the order of b. When Rc/b 
is small, we should recover the behavior in absence of ob- 
stacles. If on the other hand Rc/b is sufficiently large, 
a new regime dominated by the obstacles should appear. 
In Fig. [TJ we present our data for (R 2 )/b 2 as a func- 
tion of the scaling variable x = (N)/b 1 ^ f . The results 
refer to the trefoil knot (n c = 3). The data for various b- 
values collapse on a single curve, which can be fitted as a 
combination of two straight lines, the slopes of which are 
respectively 2v = 1.15 ± 0.02 and 2v = 0.99 ± 0.02. The 
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FIG. 1: (Color online) Log-log plot of the mean squared radius 
of gyration scaled by b 2 as a function of x = N/b 1 ^ ' 589 for 
three different b- values (results for trefoil knot). The dashed 
line has a slope 1.16, the dot-dashed one 1.00. 

former value equals (within the accuracy) the expected 
value 2vf. The latter coincides precisely with that for 



a three dimensional branched polymer (BP) for which 
vbp = 1/2 [23]. That unknotted ring polymers in an ar- 
ray of obstacles should behave as branched polymers was 
predicted in Q. Our results show that the same is true 
for knotted polymers and provide the first description of 
a polymer that is both branched and knotted (Fig. 
The crossover that we find here was never observed be- 
fore. It occurs around x = x c ~ 20. Fig. Q] clearly shows 
that at the crossover, Rg ~ b. 



FIG. 2: (Color online) Branched polygon with a composite 
3i(| 3i knot in an array with 6 = 7. 

Another important property of knotted polymers is the 
size 1^ of the knotted region. Intuitively it corresponds to 
the minimal number of monomers that still contain the 
topological entanglement of the given knot type (for op- 



erational definitions of see [18J). Research in the last 
decade has clearly shown that the size of a knot is a 
fluctuating quantity whose statistical properties strongly 
depends on the physical regime in which the polymer ex- 
ists . In good solvent conditions, it was found that 
a knot is weakly localized. This means that the average 
number of monomers in the knot (Ik) grows as a power of 
N, (Ik) ~ N f , with an exponent t < 1 Numerically, 
one has t ~ 0.75. In contrast, in the collapsed globular 
phase, the knot delocalizes (t = 1) [18l. \26\. 

For the knotted polymer inside the array of obsta- 
cles we determined Ik using the 'cut and join'-algorithm 
presented in [l8| . In Fig. [31 we report our results for 
{Ikj/b 1 ^ 1 as a function of x for the trefoil for different 
6-values. The results collapse once more on one curve. 
Moreover, there is again a crossover between two regimes. 
In the first one, the knot is weakly localised with an ex- 
ponent that we estimate as t = 0.69 ± 0.04, consistent 
within error bars, with the value in good solvent [l8j]. At 
higher values of x, the exponent t = 1.05 ±0.15 indicates 
a full derealization of the knot. It is remarkable that 
this delocalisation crossover occurs at a higher value of 
x than that of the geometrical crossover, namely near 
x = Xd ~ 100 . This suggest that the delocalisation of 
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FIG. 3: (Color online) Log-log plot of the mean knot size 
scaled by b 1 ' "f as a function of x for three different 6-values 
(results for trefoil knot). The blue dashed line has a slope 
0.69, the purple (dot-dashed) one 1.00. 



the knot only occurs when the average squared radius of 
gyration of the knotted region, {R 2 ,) becomes of the order 
of b 2 . In Fig. @]we have therefore plotted {R 2 )/b 2 versus 
(lk)/b 1 / Uf . These results indeed show that the knotted 
region behaves in a way completely similar to the whole 
polymer: once Rc,k = {RlV^ 2 becomes equal to the dis- 
tance between the obstacles, its geometry changes from 
self-avoiding to branched. Comparing Fig. [3] and Fig. 
IH we see that the knot indeed delocalises when i?G,fe is 
close to b. 
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FIG. 4: (Color online) Log-log plot of the mean squared radius 
of gyration of the knotted region (scaled by b 2 ) as a function 
of (ifc)/& 1//0 ' 589 for three different 6-values (results for trefoil 
knot). The blue dashed line has a slope 1.18, the purple (dot- 
dashed) one 0.98. 

To summarize, there are three regimes for a knotted 
ring polymer in an array of obstacles. In the first regime 
the knot is weakly localised and geometrically both the 



knot and the whole polymer behave as in good solvent 
with v = Vf. The polymer is essentially confined in one 
single cell of the array of obstacles. In a second regime, 
for 20 < x < 100, the polymer spreads over several cells 
but the knot is still confined to a single one. As a con- 
sequence the geometry of the whole polymer becomes 
branched, whereas the knotted part still behaves as in 
good solvent. In these two regimes, the knot is localised. 
Finally, for x > 100, the knot delocalises and assumes a 
branched shape just as the whole polymer. 

We next investigated the dependence of our results 
on knot complexity by extending our calculations to the 
prime knots 4i and 6i and to the composite knot 3 if) 3i 
and, for comparison, to the unknot. For this study we 
fixed 6 = 7. Firstly, we found that K c did not depend 
on the type of knot considered. This result was already 
established for knotted polygons in free space and our 
result shows that it is also true within a gel. 

In Fig. [5] we present our results for (R 2 )/(N) as a 
function of (N) for the various knot types. Plotted in 
this way, the curves reach a constant in the BP-regime. 
This figure shows several interesting features. Firstly, one 
observes that for small (N), the average squared radius 
of gyration depends on knot type. The more complex the 
knot, the smaller is the polymer. This is a result similar 
to that found for ideal knots [27|- Secondly, since more 
complex knots are smaller they feel the presence of the 
obstacles only for larger (N) . Hence, x c depends on knot 
complexity. Finally, one notices that in the BP-regime, 
polymer size hardly depends on knot complexity. This 
result is rather unexpected. Therefore in the BP-regime 
one can hardly determine the knot topology from geom- 
etry. This in contrast to the good solvent regime where 
the part of the polymer occupied by the knot, as mea- 
sured by {Ik), increases with knot complexity. Evidence 
for this can be seen in Fig. [6] where we plot our results 
for {lk)/lmin versus (JV). Here, l m - ln is the minimal num- 
ber of edges necessary to embed a given knot type in the 
cubic lattice. For the trefoil, the 4i and 6i-knot, l m i n 
equals respectively 24,30 and 60 [28j. The data shown 
in Fig. [5] show that more complex knots occupy a larger 
fraction of the polymer. Remarkably, there seems to be 
an exact proportionality with the minimal length. 

Our results can have implications for gel electrophore- 
sis experiments. Assuming as a first approximation that 
the knotted ring polymer is spherical with radius Rg, 
the mobility in an electric field /i ~ Rq 1 , should depend 
on knot type for small x- values. In this regime, different 
knot types can be separated by electrophoresis. How- 
ever, our results also predict that for x sufficiently large 
this is no longer possible. The genome of the P4-phage 
is an 11.6 kbp DNA. Its size has been estimated from 
atomic force microscopy images for DNA with complex 
knots and was found to be about 300 nm (29[, which 
is of the same order as that of the holes in an agarose 
gel. Similar considerations can be made for the knotted 
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FIG. 5: (Color online) Average squared radius of gyration per 
monomer, (R 2 ) / (N) versus average number of monomers for 
different knot types in an obstacle array with 6 = 7. 



10000 



1000 



V 




1000 
<N> 



10000 



le+05 



FIG. 6: (Color online) Average knot size divided by minimum 
size for the given knot type as a function of {N) for three 
different knot types (b = 7). The two straight lines have 
respectively slope 0.69 (orange dashed) and 1.0 (magenta dot- 
dashed). 

7 kbp DNA of plasmids. While x c is of course not a 
universal quantity, our rough estimates do suggest that 
the behavior found here is relevant for electrophoresis 
experiments. Microlithographic arrays of self-assembled 
magnetic beads pjj could also be used to reveal the 
crossovers. 

To summarize: for knotted polymers in an array of 
obstacles, we have found three regimes with distinct ge- 
ometrical and topological properties. These regimes are 
manifestations of different degrees of localization of the 
topological entanglement within the polymer and the ar- 
ray of obstacles. 
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